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ABSTRACT 


This investigation deals with distortional buckling of 
non-pri smat ic cantilever beams. A finite element formulation is 
developed and its convergence and accuracy are demonstrated. Elastic 
buckling solutions are obtained for a mono-symmetric cantilever beam 
under a concentrated load at its end. The effect of both the flange and 
the web tapering along with the effect of the application of load at 
top flange/ shear centre on the buckling load are studied. The 
influence of different types of braces and their location is examined 
in detail on distortional buckling of beam for various degrees of 


mono-symmetry . 
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CHAPTER I 

INTRODUCTION 


1.1 Stability 

The stability of structure essentially means the stability of 
its equilibrium configuration or state. In mathematical sense, 
stability is usually interpreted to mean that the infinitesimal 

disturbance will cause only infinitesimal departure from the given 
equilibrium position. In practical sense, an equilibrium state of a 
structure or a system is said to be stable if accidental forces, 

shocks, vibrations, eccentricities, imperfections or other possible 

irregularities do not cause the system to depart excessively from its 
designed position. 

In designing structures, care should be taken that the 
stresses developed do not exceed certain limits, which may otherwise 
lead to failure. Two kinds of failure can be associated with the 
structure (i) failure with respect to material behaviour, (ii),form 
failure. In the first case, stresses exceed the given safe limits 

'resulting in formations of cracks which cause failure. In the second 
case, the stresses do not exceed the safe values but structure may not 
keep up the forms in which it is designed. Also, loads on a structure 
could be tensile and/or compressive. The loss of stability due to 
tensile loads falls in the broad category of material instability, 
whereas the stability loss under compressive load is termed as 
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structural ( or geometric or form ) instability, commonly known as 
buckl ing . 



[a) Local Buckling 



Fig. 1.1 MEMBER BUCKLING MODES 

Consider the case of an I-beam under the action of transverse 
load. The beam may buckle in three ways , ( i ) local buckl ing, ( i i ) 
lateral buckling, (iii) distortional buckling, as shown in Fig. 1.1. 
The local buckling is characterized by the change in cross-section 
geometry without overall lateral displacements or twist. It can be 
further subdivided into flange buckling , web buckling or coupled 
buckling (in which both the web and the flange buckle together). 
Lateral buckling involves lateral displacements and/or twist without 
local changes in the cross-section geometry. It is also possible for a 
member to buckle into a mode which combines lateral displacement and 
twist, together with local changes in the cross-section geometry. This 
type of buckling is as termed distortional buckling. The web 
distortion which allows the flanges to deflect laterally with 
comparatively little twist, reduces the effective torsional resistance 
of the member and consequently reduces the resistance to 
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f 1 exure-tors iona 1 buckling. The effect has been found to be more 
significant in plate girders with slender, unBtiffned webs. More 
generally, distortional buckling occurs in a mode which is intermediate 
between the local buckling (which occurs in very thin webs) and the 
lateral buckling (which occurs in members of thicker web by ignoring 
any distortion). 

When a thin walled member of mono-syunetr ic I-section 
twists, the longitudinal normal stress due to the axial force and/ or 
the moment exerts torque about the axis of twist. This disturbing 
torque is commonly referred as the Wagner Effect, which results in an 
effective increase or decrease in the torsional stiffness depending on 
whether the axial force is tensile or compressive, and whether the 
smaller flange is in compression or in tension. 

The elastic buckling capacity of a member may be greatly 
enhanced by the introduction of some forms . This could be achieved by 
one or more discrete braces, or by a continuous medium such as a 
diaphragm. A discrete brace may act to restrain either the lateral 
displacement or the twist or both in the beam cross-section at the 
brace point. It may also provide some minor-axis rotational or warping 
restraints. Most of the bracing arrangements IFig. 1.2 a] in an I-beam 
may be represented by an idealized spring system consisting of an 
elastic lateral brace at a certain distance above the shear centre of 
cross-section and an elastic rotational brace IFig. 1.2 b] . A full 
brace, as oppose to a partial brace, prevents beam from lateral 
translation and rotation in the cross-sectional plane at the braced 


point . 
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(a) Typical Braces 



(b) Idealized Braces 


Fig. 1.2 TYPES OF BRACES 

1.2 Techniques For Buckling Analysis 

There are a large number of techniques available which can be 
grouped under two general approaches viz. the static-kinematic approach 
and the work-energy approach. The two approaches correspond to the 
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different strategies used in satisfying the state of neutral 
equilibrium for the deformed member. 

In the static-kinematic approach, the field differential 
equations are derived from (a) the assumed stress-strain laws, (b) 
equilibrium equations and (c)strain-displacement relations . The last 
two conditions are obtained from static and geometric considerations on 
a free body of the deformed member. Only in some rare cases closed 
form analytical solutions are possible, whereas using numerical 
technique such as the finite difference method, the finite integral 
method and the Galerkin’s method one can obtain reasonably accurate 
resul ts . 


In the work energy approach the equilibrium condition of the 
deformed member is satisfied via the virtual work concept. This can be 
done by establishing the potential energy functional and assuming the 
external forces to be unchanged during displacement variation. The 
important principle which applies for elastic problems is the principle 
of energy conservation. On the basis of the principle of energy 
conservation and the principle of stationary potential energy, various 
energy methods were proposed by early researchers e.g. the 
Timoshenko-energy method, the Ray leigh-Ritz method and the finite 
element method (FEM) . The FEM is an extension of the Reyleigh-Ritz 
method in that it involves the division of the structural member into 


member elements. 



1.3 Literature Review 


/ 


The early work on elastic buckling of I-beams remained 
confined to doubly-symmetr ic cross-sections. This led to several 
theoretical and experimental studies resulting in the development of 
approximate formulae and design charts which are now available in 
standard texts, e.g. Timoshenko (1961), Trahair (1988). Anderson and 
Trahair (1972) presented results for simply supported and cantilever 
mono-symmetric I-beams with concentrated and distributed loads and 
investigated the influence of the load height above the shear center on 
the elastic buckling moment. They proposed some modification factor to 
approximate the maximum moment in terms of the critical moment for 
buckling of a doubly symmetric I-beam under a moment gradient. This, 
and similar other studies on mono-symmetric beams, highlighted a 
difficulty : how best to determine the mono symmetric section property 
so that the flexure-torsional buckling moment could be expressed in a 
convenient form. This difficulty was overcome by Lit iporancha i and 
Trahair (1980) by introducing a new mono-symmetric parameter and thus 
simplifying the complicated expression being used in earlier 
investigations. This work led to more detailed studies of 
mono-symmetric beams by Vang and ki t iporanchai (1986) for a simply 
supported beam under a uniform moment. Subsequently Wang and 
Kit iporanchai (1989) introduced few more additional parameters and 
investigated the mono-symmetric effect under various loading 
conditions . 

For tapered mono-symmetric I— beams, Kitiporanchai and Trahair 
(1975) showed that the shear centre axis is not parallel to the 
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centroidal axis. Thus the minor axis bending moment and the torsion can 
not be separated. A finite element formulation for such beams was 
developed by Bradford (1988). He investigated the effect of flange 
tapering and web tapering on the elastic critical buckling loads 
for various degree of mono-symmetry. Kumar and Shah (1991,1993) also 
studied the effect of two way tapering of simply supported I-beam and 
beam-columns for various loading conditions and compared the results 
with those obtained using IS:800 specifications. 

The first study on distortional buckling of cantilevers in 
the open literature is available due to Johnson and Will (1974). In 
this study two relatively fine meshes of plane-stress bending finite 
elements were used. Rectangular plate elements were used to model a 
cantilever beam loaded at its end and a simply supported beam loaded at 
midspan by a concentrated load applied at the top flange. These results 
were later verified by Hancock and Trahair (1978) using energy method. 
They also presented some approximate closed form solutions which could 
be used for design purposes. Hancock (1978) investigated the buckling 
of an I-beam bent about the major axis and determined the interaction 
between the local distortional and the flexure-torsional modes. Later 
Hancock and Trahair (1980) used the finite strip method to study the 
effect of web-distortion oh elastic-flexure torsional buckling of 
simply supported beam-columns in uniform bending; a closed form 
solution was obtained as an interaction equation between the critical 
moment and loads. Bradford (1990) studied the effect of distortion on 
beams of tee-section under a uniform moment using a fine mesh of 
elements for discretization of beam. Bradford and Trahair (1981) 
developed a line finite element method for the analysis of distortional 



buckling of doubly symmetric I-beams. This method used line elements 
for beam discretization but allowed the cross-section to distort. 
Later, this formulation was modified by Bradford (1992) to take into 
account the effect of the load height above shear centre. 

Among the early studies on elastic buckling of beams with 
discrete braces, the study by Flint et al (1973) presented theoretical 
stiffness requirement for beams and columns with midspan lateral 
braces. Much later, Ki t iporanchai et al (1984) investigated the effect 
of providing braces on the critical load of cantilever beams. But, the 
most notable reference is that of Wang et al (1986) on elastic lateral 
buckling of mono-symmetric cantilevers. In this investigation the 
influence of lateral, rotational and the full brace at various 
locations along the length of beam and of the load-height above shear 
centre was demonstrated on the critical buckling moment. Bradford 
(1988) developed the finite strip method to study the effect of web 
distortion on the critical load of mono-symmetric I-beams with 

stiffened webs. However, this method proved to be cumbersome for such 

\ 

appl icat ions . 

1.4 Scope And Aim Of Investigation 

Modern methods of fabricating steel members allow welded 
I-beams to be made easily, and it is often economical to produce such 
beams with thin unstiffned webs and unequal flanges. In such cases, the 
distortion of the cross-section is likely to modify the elastic 
buckling stress of beams. The elastic buckling stresses are important 



for design since these are directly used in obtaining maximum 
permissible stresses. 

As has been pointed out in the previous section, most of the 
studies concerning cantilever beams have been confined to invest igat ing 
the 1 ateral -tors ional buckling of symmetric/mono-symmetric I-beams, 
with or without web/flange tapering. Investigations on distortional 
buckling are rather few and these too have been restricted to symmetric 
I-sections. Except the work of Bradford (1980) no study appears to be 
available on the distortional buckling of mono-symmetric sections. 

While designing a cantilever beam it could be economical to 
design it by taking into consideration the tapering of its web and/or 
flange. Moreover, the buckling capacity may get enhanced by the 
introduction of some form of intermediate braces. Hence, this 
investigation aims at studying the distortional buckling of mono- 
symmetric, braced, cantilever b6ams with web and flange tapering. It 
makes an attempts to 

* study the effect of the mono-symmetry on the buckling load of 
different cross-sections for various beam parameters, 

* study the effect of web and flange tapering on distortional 
buckling of I-beams, 

* study the effect of various braces (i.e. the lateral brace, tbe 
rotational brace, the full brace) and their location on bucking 
load of cantilevers, 

* to demonstrate the significance of incorporating distortional 
buckling as against studying lateral buckling only by treating 


the section as rigid. 



CHAPTER II 


PROBLEM FORMULATION 


2 . 1 General 

In this chapter, the finite element formulation is presented 
for analyzing the distortional buckling of mono-symmetric cantilever 
I-beams. The method is similar to that described by Bradford and 
Trahair (1981) and by Bradford (1991) for prismatic members having 
doubly symmetric sections. But, in using uniform elements to model a 
tapering mono-symmetric beam, one may expect a slow rate of 


convergence . 

It 

is. 

therefore , 

desirable 

to develop 

a tapered 

finite element 

that 

accurately 

represents 

a tapered 

member and 


does not introduce any artificial discontinuities at node points. 
The method requires a lengthwise subdivision of beam into one 
dimensional elements, for which the stiffness equations for buckling 
are developed. 

The axis system of each element is shown in Fig. 2.1(a). The 
z-axis is the longitudinal axis, while the x- and the y-axes are the 
major and the minor principal axes, respectively. The actions of 
various loads on the beam element are shown in Fig 2.1(b), and are 
applied initially in the plane of symmetry y-z . The initial action set 
is W1 , Ml and VI at the end 1 (i.e. z=0) and W2, M2 and V2 at the end 
2 (i.e. z=L); these actions are increased proportionally by the load 
factor X until buckling occurs. 




2.2 Section Properties 


Following parameters are found convenient to define the section 
properties : 


(i) Mono-symmetry parameter (p) 


P 




I 


yb 


( 2 . 1 ) 


where I 4 and I are the moments of inertia about the y-axis of the 

y t yb 

top and the bottom flanges, respectively. The parameter p=l represents 
a T-section, p = 0.5 represents a doubly symmetric I-section, whereas 
p=0 corresponds to an inverted T-section. 


(ii) Polar moment of inertia (J) 

for thin open cross-sections it can be approximated as 


J 


n 

; £*, 


i * 1 


( 2 . 2 ) 


where b ^ and t. are the length and the thickness of i th member, 
respect ively . 

(iii) Beam parameter ( k ) 



„2 2 
n e i h 

y 


4 G J 1 


(2.3) 


where I is the moment of inertia of the whole section about the 
y 

y-axis and 1 is the length of the beam. E and G are the modulus of 
elasticity and the modulus of rigidity, respectively. 
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2.3 Longitudinal Variation of Displacement 


The cross-sect ion 

of the beam i s 

assumed to 

cons i st 

of 

two 

rigid flanges connected by 

a flexible web 

which may 

d i stort 

in 

the 


plane of the cross-section. The displacement of a cross-section of a 
beam element is shown in Fig, 2.2 . Vhen the beam buckles, its flanges 
deflect u and u. rotate u and u , and twist <P and 4> , where a 

t b t b t b 

prime ( )' indicates differentiation with respect to z. The lateral 
displacements of the flanges are assumed to be cubic polynomials, so. 
that 


u t 

= a j L + 

a 2 L 

z 

+ a 3 L 

z 2 

+ « 4 L 

7' 

u b 

= 0£ L + 

5 

*6 L 

Z 

+ “7 L 

z 2 

+ a e L 

z 3 


where Z = z/ L, and L being the length of elements. The flange twist 
may be expressed as 


4> 

t 


u 


t 


u 


b 


h 


+ a 


9 



4 > 

b 



+ a + 
1 1 


a 


\ z 


Z 


(2.5) 


in which h is height of the section. Eqs . (2.4 ) and (2.5) may be 
expressed in matrix form by 


{ u ) 


I M ] .{ « ) 


( 2 . 6 ) 
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in which 


{ U } = { U t , U b , H> t , Q h 
{ “ > = { a ,- tt 2 '“a 


T 

) 


a io' a t 1 


<S 2 } 


and [M3 is given in Appendix -I. The nodal displacement vector 
can be represented as 


{ 5 } 



■ ■ ■ V 

U t2’ U bl' U b2’ U bl' U b2'* > tl'^t2'* > bl’^ > b2j 


which represents the degree of freedom of nodes 1 and 2 
element. By using eq . (2.7) the vector { 6 }can be expressed as 
{ d ) = [ C ) { a } 

and so 

{ « } = l C j" 1 { 6 } 

The matrix [ C ] 1 is given in Appendix -II. 


2.4 Cross-sectional Variation of Displacements 


Since the flanges are assumed to be rigid, the vertical 
displacement of any flange point x is either 

v = x <t> L or v = x 

it b b 


(2.7) 

$ } 

( 2 . 8 ) 

of the 

(2.9) 

( 2 . 10 ) 


( 2 . 11 ) 


Moreover, it is assumed that the web becomes a cubic curve; thus 
the deflection at any point y can be written as 


( 2 . 12 ) 


u * = < “,3 + «14 Y + “« Y “ + a t6 Y " J f ( z > 


in which Y = 2 y/ h and f is a cubic function of z . Using the 
ident it ies 


u « ( u ) 

t v w -h/2 


u b = ( u h/2 


V 


V 


9 u 


9 y 
9 u 


d y 


J -h/2 


h/2 


C2. 13) 


the web displacement can be expressed in terms of the vector of 
nodal displacement { 5 ) by 


-l 


u = {N} [Cl {6} 


(2.14) 


The vector 


-a 


{ N } = { N } 1C] [ M ] 


(2.15) 


is given in Appendix -III 


2.5 Element Stiffness Matrix 

During buckling, the beam element deflects and twists. The 
increase U in the strain energy during buckling can be approximated by 

U = U„ + U + U (2.16) 

f U f^P W 

9 

in which 

U « the flange in-plane bending strain energy. 


r ei 

f u" ] 

+ EI _ 

fu" ] 

J yft 

l t J 

yfb ' 

l b J 


1 

2 


dz 


(2.17) 



1 


where I and I yfb are the moments of inertia about the y-axis of 
top and bottom flanges, respectively. 


U t( p = the flange out of plane bending and twisting energy 


* 5 D , 


b / 2 


[ 1 1 ( v t" )*+ 2 <» - « > ( y**} 

&£ QZ 


-b t /2 


b /2 
b 


. J ”( v b' )"♦ 2 <> - « > ( ~z~r~ ) 2<1 * ] <2 - 


•b /2 

V 


&z 


18) 


3 Z 

where v is the Poisson’s ratio and D f - E T / 12 (l-x> ) and b t and 
b L are the width of top and bottom flanges, respectively. 

b 


Uw = the 

strain energy of 

web 




r L h / 2 

r 1 



x 1 (\J 

11 


J < e 3 

{ o } dy dz 

(2.19) 


o J 

- h/2 




In eq. (2.19), the web generalized strain vector { e } is 


r 2 

/ a u „ 

z 

l U W - 

b\ i 

(2.20) 

X 2 

d y 

w 2 ’ " 

b Z 

by bz ' 



while the generalised stress vector is given by 

i 

{ o ) = ID 1 { e ) (2.21) 

w 


in which 



1 


V 


0 


I D 


w 


E 

12 


(1-V) 


v 1 

0 0 


0 

(1-u) 

2 


( 2 . 22 ) 


The stiffness matrix I k ] of a beam element can be derived by 
expressing the strain energy in the form 

V = l { 5 } T [ k ] { 6 } (2.23) 

Substituting eqs.(2.17) into eq . (2.16) and noting that eq . (2.16) and 
eq . (2.23) apply for all { 5 }, lead to stiffness matrix 

l k ] = [ C )“ T I k a ] [ C l" 1 (2.24) 

The procedure of calculating the elements of stiffness matrix I k ] Is 
explained through an example in Appendix -IV. 


2.6 Element Stability Matrix 

The beam is subjected to bending moment and shear forces. 
The decrease in the potential energy of these actions during buckling 
is given by 


V = V 


f u 


+ V 


f«> 


+ V 


(2.25) 


in which 


V, « the work associated with the in-plane flange 
fu 

deformation 



{ u* }' 


fb 


{ U* } 

O 



dz 


1 

2 


(2.26) 



where A ft , A fb are the areas of and o^ t are the bending 

stresses in the top and the bottom flanges, respectively, 


V r4> = the work done associated with the out-of plane 
flange deformations 

l f L b t /2 b b /2 

' 5 T J [ J V v i ,Z * J % < v : 3 2 * 1 ] <lz (2-27) 

-b t /2 -b b /2 


V • the work associated with the web deformations 
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? dy dz (2.28) 


M 


(here o = X 


y , 


T = x 


h.t 


(2.29) 


In eq.(2.29) M is the moment about x-axis, V is the shear force at 
x 

that section and I is the moment of inertia about the x-axis. 

x 

The element stability matrix X I g ] is obtained from the 
decrease in potential energy V by a procedure analogous to the one 
described for the stiffness matrix in the section 2.6, and is given by 


X [ g ] = X [ C ]“ T I 1 IC1 


(2.30) 


2.7 The Work Done by Loads not at Shear Centre 


It is well known that when a concentrated load V acts 



at a height a above the shear center of the beam, the lateral 

buckling resistance is lowered because of the additional torque W,a.<J> 

which increases twisting of the beam. When the load XV at element 

1 

node 1 acts within the web height a above the shear center the 
point is lowered by ^ 0 J a ( fcu^ / &y ) 2 ; hence the loss of potential 

energy is given by 



-h/2 < h/2 (2.31) 


The inclusion of the load height in finite element model was achieved 
by calculating and adding separately these terms to stability matrix 
X [ g J. 


2.8 Effect of Mono- symmetry 

When a beam of mono-symmetric section is loaded in its 
plane of symmetry, it twists during buckling and the longitudinal 
bending stress M^y /i ^exerts torque. The action of this torque is 
to change the effective rigidity of the beam. The compressive bending 
stress causes a disturbing torque while the tensile stress a gives a 
restoring torque. In a doubly symmetric section, these two torques 
balance each other. On the other hand, in a mono-symmetric section 
there is an imbalance torque which depends upon the lever arm between 
the flanges and the shear center. Since in the present investigation 
the web distortion is being taken into consideration, the location of 
the shear centre is a difficult proposition. Hence, for calculating 
the work done by the imbalance torque, it is assumed that the shear 
center lies on the straight line joining the top and the bottom 
flanges. Accordingly , this work done can be expressed as 
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X M P ( <t> ) dy 

x x a v g 


( 2 . 32 ) 


in which 


„ l f 2 3 ,, 

**x = 77- J ( * y + y ) dA + 2 y , 


4> = ( u - u .) / h , 

( a v g . ) i b 9 * 


and y is the depth of shear center from the top flange. 


2.10 


Global Stiffness and Stability Matrices 


The global stiffness matrix [ K ] and stability matrix [ G ] 
can be assembled from the individual element stiffness and stability 
matrices [ k ] and [ g ], respectively, using equilibrium and 
compatibility at nodes. The numerical integration with respect to x.y 
and z is carried out using the four-point Guassian-quadrature . The 
total potential 0 may be written as 


n = | t Q ) T [ I K ] - X [ G ] ] { Q ) 


( 2 . 33 ) 


where { Q } is the vector of global degrees of freedom assembled 
from the { q } vectors of each element. Invoking the principle 
of stationary potential energy yields 

( 2 . 34 ) 


511 = { SQ ) [ IK] - X [G] J {Q> = 0 


Noting that the variations { 5 q } are arbitrary, this results in 


[lKj-XIG)]{Q}=0 1 2 - 35) 

Mo. h. J 



2 ; 


The values of X which yield a non trivial solutions for { Q } in 
eq. (2.35) are the eigen-values, while the corresponding values of 
{ Q } are the eigen-vectors. For a non-trivial solution 

| [ K 1 - X t G ] =0 (2.36) 

The above eigen-value equation is solved using the ’NAG-LIBRARY’ 
subroutine *F02BJFE.f’. The critical load is given by lowest 
positive eigen-value. 
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CHAPTER III 

RESULTS AND DISCUSSION 


3.1 Introduction 

The finite element formulation developed in chapter II is 
used to study the 1 atera 1 /di stor t iona 1 buckling of mono-symmetric, 
tapered, braced cantilever I-beams, under a concentrated point load at 
the free end. As mentioned earlier, the main assumption in the 
formulation has been that the flange remains undistorted during 
buckling. Because of this assumption the local buckling for low h/t 
ratio ( web height to thickness ratio) can not be predicted accurately 
( Hancock, 1978,1980). 

Since tabulated closeddorm solution for distortional buckling 
are not available for mono-symmetric, tapering or braced cantilevers, 
the stiffness and the stability matrices generated are tested by 
comparing value of buckling load with available results ( Bradford, 


1988; Trahair, 1986) 

for beams with rigid 

web . 

To 

reduce the 

distortional buckling 

case to 

the rigid web 

case 

C 

i.e. lateral 

torsional buckling ) 

the plate 

rigidity, , 

for 

the 

out-of-flange 


flexure of web is made to approach infinity. 

The geometric parameters of the beam cross-section are shown 
in Fig. 2.1. In presenting the results of the effect of mono-symmetry 
on the distortional buckling, the distortional load, is normalized 
with respect to the lateral buckling load, W or ,( i.e. considering the 



24 


web as rigid). On the other hand, in discussing the effect of flange 
and web tapering on distortional buckling, the load is 

non-d i mensi onal i sed W 1 2 / ■/ (El GJ) o in terms of beam parameters. 

Throughout, the value of Young’s modulus is taken as 2 X 10 8 kN/ m 2 and 
that of Poisson’s ratio as 0.3. 

3.2 Accuracy and Convergence 

The distortional buckling of cantilever beams with 
concentrated load at the shear centre of the free end was considered by 
Johnson and Will (1974). The geometrical parameters were : h = 533 mm., 
b fa = b t = 254 mm, 1 = 3810 mm. In their analysis, two relatively fine 
meshes of rectangular elements were used. One mesh consisted of 60 
elements (20 in the web, 20 in each flange), while the second consisted 
of 96 elements (48 in the web and 24 in each flange). The buckling 
loads for the two cases were 1410 kN and 1370 kN, respectively. Using 
the present formulation, the buckling loads are 1413 kN ,1405 kN and 
1403 kN for 4,8 and 10 longitudinal elements, respectively. The lowest 
value (1403 kN) differs by 2.4 % from that of Johnson and Will (1974). 
The rate of convergence of the finite element method is investigated 
for a doubly-symmetric beam, both with tapered web ( 0 = 0.5) and with 
tapered flange( a = 0.5 ), with point load applied at the shear centre 
and at the top flange. 

Fig 3.1 shows the convergence of the non-dimensional load, 

W l Z //(EIGJ)o, with the number of elements ( the subscript o refers to 
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Shear Centre Loading Top Flange Loading 

Fig. 3.1 CONVERGENCE OF BUCKLING LOAD FOR CANTILEVER 
( ) Tapered Element , ( ) Stepped Element 

the values at the fixed end ). It can be seen that the convergence is 
very rapid in case of tapered element; the load remains more or less 
constant for four or more elements. Moreover, the convergence is faster 
while considering the effect of the flange taper than that of the web 
taper. It is also observed that the top flange loading requires 
comparatively more number of elements as compared to the shear centre 
loading. In any case, minimum of ten elements are used to model the 


beam . 
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3.3 Effect of Mono-symmetry 




b fc /T = 10 b fa /h =0.2 b fa /T = 15 b^h =0.4 


Fig 3.2 EFFECT OF MONO- SYMMETRY ON DISTORTIONAL BUCKLING. 

(Shear Centre Loading) 

( ) h/t =1 00 ( ) h..'t=2S0 

Figures. 3.2 and 3.3 show the effect of mono-symmetry on the 
distortional buckling of I-beams. Similar to the case of 

lateral-torsional buckling, the beam having the bigger flange in 
compression giveB higher buckling load than the beam with equal flanges 
or with smaller flange in compression, for the same value of the beam 
parameter K. To demonstrate how the distortion affects the buckling 
load when flanges are unequal, the distortional buckling load , ^ od > is 




27 


normalized with respect to latera 1 -tors ional buckling load, V , and 

or 

the results are plotted against 1 /h ratio for different values of the 
mono-symmetry parameter p. Since the compression flange width to the 
height ratio ( b / h ) is also one of the influencing parameter, the 

D 

tension flange width b varied (keeping the compression flange width b 

1 b 

constant) in order to get various values of p. 




b b /T = 10 b b /h =0.2 


b fc /T = 15 b b /h =0.4 


Fig 3.3 EFFECT OF MONO-SYMMETRY ON DISTORTIONAL BUCKLING. 

(Top Flange Loading) 

( ) . h..''t-=100 , ( ) h/t =250 
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Figures 3.2 and 3,3 show that as the degree of mono-symmetry 
decreases, the reduction in the buckling load ratio also decreases, 
i.e. bigger the compression flange, lesser is the distortional buckling 
load. Further, as the b /h ratio increases, the effect of distortion on 

D 

buckling load becomes prominent. It can also be observed that the 
effect of web distortion becomes pronounced with the increase in the 
h/t ratio i.e. when the web becomes more slender. Moreover, as the 1 /h 
ratio decreases, the distortional buckling load gets reduced 
signi f icant 1 y . 

For the top flange loading, although the nature of the curves 
are similar ( Fig. 3.3 ) the effect of mono-symmetry is more than that 
in the shear centre loading. Also, the trend in the result as the 
degree of mono-symmetry varies from p =0.1 to 0.9 is quite different 
than what is observed when load is at the shear centre. 

3.4 Effect of Flange / Veb Tapering 

Two cases of tapering are considered : <i) flange tapering, 
measured in term of parameter « C ratio of the flange width at the free 
end to that at the fixed end), (ii> web tapering, measured in term of 
parameter P (ratio of the web height at the free end to that of at the 
fixed end). Fig. 3. 4-3. 7 illustrate the effect of web and the flange 
tapering on distortional buckling for the shear centre loading and the 
top flange loading. The non-dimensional ised buckling load 
W l 2 / / (EIGJ)o, is plotted against the beam parameter 

K f« D / (EI y h 2 /4GJL 2 )oj . For the sake of comparison, the corresponding 
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resultB for lateral buckling are also given in all figures. 



h/t=100 h/t=250 

Fig. 3.4 EFFECT OF FLANGE TAPERING ON DISTORTIONAL BUCKLING 

( Shear Centre Loading ) 

C ) Lateral-Torsional Buckling, ( ) Distortional buckling 

Fig. 3.4 shows the the effect of the flange tapering on the 
distortional buckling load for the shear centre loading. As the flange 
tapering increases ( « changes from 1.0 to 0.25 ) the net decrease in 
the load also goes down i.e more the flange tapering, less will be the 
distortional effect. The h/t ratio also affects the distortional 


buckling capacity for given values of flange tapering («) of beams; the 
loss in strength is more for higher h/t ratios. 




W L Z //<EIGJ) 
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Similarly 

Fig. 3.5 shows the 

effect of web tapering. 

For K 

lying in 

the range 

1. 0-4.0, there is a 

significant difference 

in 

the 

loss of 

strength due to distortion. 

but for K > 4.0 the loss 

in 

strength 

for each 

value of the web tapering parameter 

0 

i s 

approximately the 

same. Further, the 

increase in the h/t 

ratio 


decreases the net buckling load capacity. 




Fig. 3.5 EFFECT OF WEB TAPERING ON DI STORT IONAL BUCKLING 

( Shear Centre Loading) 

( ) Lateral-Torsional Buckling, ( ) Distortional buckling 


From Figures 3.6 it can be seen that for the top flange 
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loading, the loss of strength due to distortion is approximately the 
same for all values of a (flange taper ing) .when the load is at the top 
flange. Moreover, the loss of strength due to distortion is more in a 
slender web ( high h/t ratio ) than in a stocky web. 




Fig. 3.6 EFFECT OF FLANGE TAPERING ON DISTORTIONAL BUCKLING 

( Top Flange Loading ) 

( ) Lateral-Torsional Buckling, ( ) Distortional buckling 


The effect of web tapering on the buckling of cantilevers 
when the point load is at the top flange can be seen in Fig. 3.7. Both 
the values of h/t, 100 and 250, are considered. The effect of 
distortion increases as the web tapering increases ( i.e. 0 decreases). 
The effect is more pronounced for higher h/t ratios (slender webs). 



W L 2 //(EIGJ> 
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A 




Fig. 3.7 EFFECT OF WEB TAPEBING ON DISTOBTIONAL BUCKLING 

( Top Flange Loading > 

— ) Lateral-Torsional Buckling, ( ) Distortional buckling 


3.5 Braced Cantilevers 

The elastic buckling capacity of mono-symmetric cantilevers 
with braces is influenced by following factors : (i) type of brace- i.e 
full brace(fb), lateral brace at the top flange (It) , lateral brace at 
the bottom flange(lb), rotational brace at the top flange(rt), 
rotational brace at the bottom flange (rb), (ii) location of brace, 
(iii) degree of mono-symmetry (v) web slenderness ratio. Due to large 
number of influencing parameters, it is difficult to present results 
comprehensively covering all the ranges. However, an attempt is made to 
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examine their effect separately. 

The effect of a particular brace depends upon how stiff the 
brace is. Normally braces having stiffness two or three times the beam 
stiffness in that direction are sufficient to reduce the displacement 
at that point nearly to zero. In this investigation, the brace 
stiffness are assumed to be infinite, which means that the 
corresponding displacement at that point is completely restrained. Five 
values of mono-symmetry parameter p are taken for invest igat ing the 
brace effect with the point load at the shear centre and at the top 
•f 1 ange . 


First, a symmetric I-section with two web slenderness ratios 
is considered to demonstrate the effect of each brace. The full brace 
appears to be more effective for all portions of the beam compared to 
other braces as it prevents all the degrees of displacement at the 
point. The lateral brace at the top flange is more effective than at 
the bottom flange. While the optimum location for the top flange 
lateral brace is near the free end, it is somewhere near the mid-span 
for the bottom flange lateral brace. Rotational braces are less 
effective than lateral braces. The top flange rotational brace is more 
effective than bottom flange rotational brace near mid-span of the beam 
while the situation is just the reverse near the free end. Further, as 
expected, the effectiveness of braces increases with the decrease in 
web slenderness ratio, i.e. efficiency is more in stocky webs than 


slender webs. 



W L Z //iElGJy 


34 



h/t=100 h/t=250 

Fig 3.8 EFFECT OF BRACE LOCATION ON BUCKLING LOAD FOR DIFFERENT TYPES 

OF BRACES ( Doubly Symmetric Section) 

Figures 3.9 to 3.13 demonstrate the effect of the load 
position when different types of braces are considered for beams having 
varying degree of mono-symmetry. As expected, the buckling capacity 
increases as the load moves towards the shear centre. This is so because 
the distortional torque reduces as the load shifts from the top flange 
towards- the shear centre. 

As is evident from Fig. 3.8 to 3.13, in general the influence 
of a discrete brace is more beneficial in cantilever beams whenever the 
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compression flange is bigger than the tension flange. In contrast, for 
p - 0.9 the increase in the buckling capacity for nearly all types of 
braces ( expect for lateral brace at the bottom flange for the top 
flange loading ) is very small. For all degrees of mono-symmetry (p = 
0.1 to 0.9) the optimum location of brace ( for all types of brace ) is 
approximately near the mid-span when the load is at the shear centre 
and it is near the free end for the top flange loading. 




Shear Centre Loading Top Flange Loading 

Fig 3.0 EFFECT OF BRACE LOCATION ON BUCKLING LOADS 

< Ful 1 Brace ) 

For the shear centre loading, the lateral brace at the top 
flange is more effective than the bottom flange lateral brace and the 
optimum position of the brace is near the free end. For the top flange 
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loading, the bottom flange rotational brace is more effective than at 
the top flange. For the top flange rotational brace, the optimum 
location lies in the range 0.71 to 0.91, whereas for bottom flange 
rotational brace it is near the free end. 




Fig 3.10 EFFECT OF BRACE LOCATION ON BUCKLING LOADS 
( Rotational Brace at the Bottom Flange > 



3 7 




Shear Centre Loading 

Fig 3.11 EFFECT OF BRACE LOCATION ON BUCKLING LOADS 
( Rotational Brace at the top Flange ) 



Fig 3.12 EFFECT OF BRACE LOCATION ON BUCKLING LOADS 
( Lateral Brace at the Bottoa Flange ) 
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CHAPTER IV 

CONCLUDING REMARKS 


The finite element method of analysis has been presented 
which is capable of making fairly accurate determination of 
distortional buckling for tapered, mono-symmetric and braced cantilever 
I-beams. The accuracy of the method has been shown by comparison with 
some known results, and its rapid convergence has been demonstrated by 
comparing it with a finite-element representation that uses uniform 
elements. On basis of this study, some salient observations and results 
can be summarized as under: 


► the reduction in the distortional buckling load is less for higher 
values of a, the flange tapering parameter. 

► the reduction in the distortional buckling load is more for higher 
values of £, the web tapering parameter; it is definitely so for k 
up to 4 . 

► the effect of distortion on the lateral buckling load decreases 
with the increase in the size of compression flange. 

► the loss of strength due to distortion is more in a beam with the 
top-flange loading than with the shear centre loading. 

► in general the advantage of braces is more when the larger flange 
is under compression. 

► although, the full brace is most effective; other braces are also 
reasonably effective when they restrain the top flange displacement 

► the optimum location of the brace is near the mid-span for the 
shear centre loading and near the free end for the top 
flange loading. 
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[ M ] 


and 


APPENDIX - I 


Matrix I M } in Equation (2.7) 



0 


0 

0 


Mu 


0 

0 

Mu 

- 1 

h 

Mu 

M$> 

0 

Mu 

- 1 

h 

Mu 

0 

M$ 


| M u j = | L, LZ, LZ? LZ 3 | 

{*}-{*•■} 

T 

{ a ) * { a i ' a 2’ a !2 } 
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APPENDIX - II 



in which h is the distance between the centroid of the top and 
the bottom flanges, and L is the the length of an element. 



i—3 f 00 


APPENDIX 


III 


Vector { N ) in Equation (2.14) 
{ N „ } = { 1 , Y , Y 2 , Y 3 } 






APPENDIX - IV 


Procedure for Calculating Elements Of Stiffness Matrix 

Using the eqs (2.6) and (2.10), we get 
{ u } = l M ] [ C J _1 { 6 } 
and hence 

D n { u } = D n l M ) l C l" 1 { 8 } 

in which D n is the differential operator of n th order. Consider the first 
term on right hand side of eq . (2.17) which is strain energy 

stored due to lateral deformation of the top flange 

. p t j 

U « - El {u")dz 

rut 2 J yft t 

0 7 

From eq (4.2) we can write 

L 2 

= 5 J E 1 [ ' M t 1 1 C rl 1 5 1 ] dZ 

L 

= l | [ t 6 } T l C ]* T [ H[ J 1 E I yft [ M * ] [ C l' 1 1 6 } J dz 
- { 5 } T [ C r 7 l f [ m' ] T E I , I m/ ) dz l C ] “ 1 { 5 } 

^ qJ *■ y it 1 

where as t M ] is the corresponding shape function for { u ). 

t. t 
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On similar lines we can evaluate second part of the right hand side of 
eq (2.17), and the right hand side of eqs. (2.18), (2.19) and combining 
them to get the element stiffness matrix [ k ]. 

£L 



